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ABSTRACT 

Using the NICMOS coronagraph, we have obtained high-contrast 2.0 /im imaging polarimetry and 
1.1 /im imaging of the circumstellar disk around AB Aurigae on angular scales of 0.3-3" (40-550 AU). 
Unlike previous observations, these data resolve the disk in both total and polarized intensity, allowing 
accurate measurement of the spatial variation of polarization fraction across the disk. Using these 
observations we investigate the apparent "gap" in the disk reported by Oppenheimer et al. 2008. In 
polarized intensity, the NICMOS data closely reproduces the morphology seen by Oppenheimer et al., 
yet in total intensity we find no evidence for a gap in either our 1.1 or 2.0 /im images. We find instead 
that region has lower polarization fraction, without a significant decrease in total scattered light, 
consistent with expectations for back-scattered light on the far side of an inclined disk. Radiative 
transfer models demonstrate this explanation fits the observations. Geometrical scattering effects are 
entirely sufficient to explain the observed morphology without any need to invoke a gap or protoplanet 
at that location. 

Subject headings: stars: individual (AB Aur) — stars: pre-main sequence — circumstellar matter — 
planetary systems: protoplanetary disks — polarization 



1. INTRODUCTION 

AB Aurigae is one of the most intensively studied 
of all Herbig Ae/Be stars, on account of its proximity, 
brightness, and youth (distance d = 144 pc; visual 
magnitude V = 7.04 mag; age < 3 Myr; spectral 
type AOe) . Steadily improving observational capabilities 
have yielded increasingly detailed views of its complex 
and dusty environment, providing many insights into 
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the nature of circumstellar disks (e.g.lGradv et al.l 119991 : 
iFukaeawa e t alj |2004t iPietu et al.ll2005D . 

In particular. lOpoenheimer et alj (|2008l hereafter 
Opp08) recently presented high-angular resolution, high 
contrast imaging polarimetry of AB Aur at 1.6 /im, 
obtained with the Lyot Project coronagraph on the 
AEOS 3.6 m telescope. These observations resolved 
the disk in polarized scattered light as close as 40 
AU (0.3") to the star. In polarized light the disk 
is not axisymmetric, but instead shows an apparent 
gap or depleted region at a radius of ~ 100 AU. 
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Such a gap may be created b y dynamical perturbations 



from forming planets (e.g. | Kuchner &: Holmanl 120031 : 
IWvattJ l2005t (jang-Condelll 120091) . Intrieuindv. Opp08 
report a faint point source within the gap, which 
they conjecture could be the perturbing object — though 
they are cautious with this identification due to its 
low statistical significance (2.8 a). The formation 
mechanism(s) of massive planets at large separations 
remain highly uncertain (|Dodson- Robinson et al.l [20091 : 
iNero fc Bjorkmanl 120091 ): any prospect for observing a 
planet forming in situ around AB Aur should be pursued 
to clarify this puzzle. 

One challenge in interpreting the data from Opp08 
is that the disk is visible only in polarized intensity, 
P = \/ Q 2 + U 2 (w here Q and U are the usual 
Stokes parameters; see Tinbcrgcn 1996 for a review of 
polarization fundamentals and notation). In their total 
intensity image (Stokes I), the residual speckle halo of 
the stellar PSF completely drowns out the disk's fainter 
light. The origin of features seen only in polarized light 
is ambiguous: any observed spatial variation may be due 
either to variation in the total amount of scattered light, 
or to a change in the polarization fraction of that light. 
The polarization induced by dust scattering depends 
strongly on the scattering angle (see Figure Q]) , allowing 
disk geometry or viewing angle to cause variations in 
the observed polarization which might be mistaken for 
intrinsic substructure within the disk. 

In particular, for AB Aur the observed celestial 
position angle of the depleted region, 333 ± 2° (Opp08) 
is precisely aligned with the disk's apparent rotation axis 
as inferred from CO emission line kin ematics (330-333°; 
iCorder et aTll2005l : IFietu et al.ll2005fl . Is this alignment 
coincidental? 

The most direct way to answer this question is to 
obtain images with enough contrast to directly detect the 
disk in total intensity, and then calculate the polarization 
fraction, p = P/I. Such observations are best obtained 
with the Hubble Space Telescope, whose ACS and 
NICMOS coronagraphs both have provided sufficiently 
high contrast to precisely and accur ately measure the 
polarization of disk-scattered light (iHines et alJ 120001 : 
iSchneider & Hinesll2007l iGraham et al.ll2007l Hines and 
Schneider, in prep). 

In this paper we present new NICMOS coronagraphic 
imaging and polarimetry of AB Aur which clarifies the 
nature of the dark "gap" observed by Opp08. These data 
were obtained as part of a coronagraphic polarimetry 
survey of young stars across a range of masses and ages 
which will be rep orted more fully in future works (see 
iPerrin et al.ll2009l for a brief overview). 

2. OBSERVATIONS AND DATA REDUCTION 

Our observation and reduction strategies follow the 
standard recommendations for NICMOS coronagraphy. 
We observed AB Aur on 2007 Sep 14 and 2007 Dec 
21 as part of program HST/GO 11155, in two visits 
identical except for a 136° difference in roll angle. After 
centering AB Aur behind the NIC2 coronagraph hole, 
nine 192 s exposures were taken using the 2.0 /im POL*L 
linear polarizers, cycling between the 0, 120, and 240° 
polarizers after each exposure, followed by one 512 s 
exposure with F110W. The telescope was then shifted to 
move AB Aur a few arcseconds away from the occulting 
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Fig. 1. — The induced polarization for some possible grain 
compositions, demonstrating that polarization depends on both 
scattering angle and dust grain properties, particularly composi- 
tion and porosity. Positive values indicate typical centrosymmetric 
linear polarization, while negative formal polarization denotes 
linear polarization oriented radially. Polarization is maximized for 
porous grains, and for scattering angles slightly above 90°. These 
models assume a power-law distribution of grain sizes from 0.03- 
200 fim with slop e -3.5, typical for YSO disks. See section 5.2 of 
IPinte et al.l ( 2008) for a discussion of dust model and computation 
details. The overplotted points show the fractional polarization 
observed around AB Aur; the inferred scattering angles are not 
symmetric due to the flaring of the disk surface, estimated ~ 10° . 
Of the models shown, the 60% porosity silicates provide the best 
fit. 

spot, after which we took two dithered 4 s exposures 
in both POL0L and F110W for photometry on the 
unocculted star. 

Achieving high contrast with NICMOS requires sub- 
tracting a color-matched point spread function (PSF), 
and for polarimetric observations, the reference star must 
be unpolarized to minimize systematic biases. (The 
direct stellar light from AB Aur should be unpolarized 
or nearly s o, due to the low l ine-of-sight extinction, 
A v = 0.25; iRoberge et al.l 120011 ) It proved challenging 
to identify PSF stars which are unpolarized yet also 
sufficiently red to match AB Aur's color (H — if =0.832), 
which is redder than the Rayleigh- Jeans slope, hence 
anything so red must either have nonthermal emission 
(often highly polarized) or else be dusty and extincted 
(likewise polarizing). After some consideration we 
identified nearby M dwarfs as the best candidate PSF 
references, and therefore observed Proxima Centauri, 
GJ 273, and GJ 447 one visit each using an identical 
observing strategy as above. We supplemented our 
program with additional PSF observations drawn from 
programs HST/GO 10847 and 10852. 

Our data reduction approach follows that of Schneider 
et al. 2005. Briefly: starting with pipeline-reduced 
images from STScI, we corrected for bad pixels and 
sky/thermal background emission in all images, and then 
median-combined the three coronagraphic images for 
each POL*L filter. The F110W coronagraphy and the 
unocculted imaging were reduced similarly. 

Obtaining optimal PSF subtraction depends on accu- 
rate registration and flux-scaling. Starting from flux 
ratios derived from phometry of the unocculted stars, 
we adopted two independent strategies to optimize the 
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subtractions: (1) A manual search visually compared 
subtractions of different reference PSFs while interac- 
tively varying the registration and scaling to minimize 
the residuals. (2) An automated algorithm evaluated 
subtractions across a grid in alignment parameter space 
for each image pair. Each subtracted image was high- 
pass filtered to reject diffuse light from the disk while 
emphasizing features with angular scale comparable to 
the diffraction limit (the characteristic size of speckle 
residuals). The variance on this angular scale was 
minimized to find the best subtraction. The best 
subtractions were obtained using GJ 273 as the PSF; our 
two optimization approaches yielded subtracted datasets 
with polarization fractions differing by < 6%, which we 
adopt as our polarization uncertainty. 

Remaining artifacts in the subtracted images such as 
diffraction spikes were masked out. The images were 
rectified for geometrical distortion, rotated to a common 
orientation, and combined to produce final mosaics in 
each filter. From these images the Stokes parameters 
I ,Q, and U were d erived using the POLARIZE software 
(jHines et al.1 120001 ) which models the imperfect linear 
polarizers in NICMOS. The output polarized images were 
then smoothed by a one resolution element (~3 pixels 
= 0.22") Gaussian, but we present the undersampled 
F110W data at full resolution. 

The resulting images are shown in figure [2j The 
F110W subtractions are excellent, yielding a clean image 
with minimal PSF residuals. The POL*L images 
are more affected by instrumental residuals due to 
optimization of the NIC2 coronagraphic optics for shorter 
wavelengths, and to a better PSF template color match 
at 1.1 fim. These factors result in 2.0 fim images dappled 
with faint residual speckles, but that still clearly show the 
bright circumstellar nebulosity. 

3. RESULTS 

A complex and asymmetric nebula surrounds AB Auri- 
gae. The F110W image clearly shows the m ultiple spiral 
arms previo usly observed, f o r inst ance by IGradv et al.1 
(|1999ft and iFukagawa et~aT1 (|2004f >. The bright inner 
region of the disk extends out to ~ 1.2", and is brightest 
to the south and southeast, as was seen by Fukagawa et 
al. At 2.0 fim in total intensity (Stokes /) the overall 
surface brightness distribution is similar to that seen at 
1 /im, albeit at lower angular resolution and contrast. 

The 2.0 fim polarized intensity image (lower left 
panel of Figure ^ reproduces with high fidelity the 
polarized intensity pattern from Opp08. Both images 
show the region of lower polarized intensity at PA=333° 
between two brighter "shoulders" on either side, the 
brightness enhancement southwest of the star, and even 
the scalloped, almost-concave southern edge of the bright 
polarized region. Comparing the NICMOS and AEOS 
datasets, the NICMOS image is more sensitive and traces 
polarized light further from the star (~ 7" vs. ~ 
1.2"), while the AEOS image has slightly better angular 
resolution, and better speckle reje ction due to the 
simultaneous diffe rential technique (jKuhn et al.1 120011 ; 
iPerrin et al.1 1200 8). Yet the most significant advantage 
of our NICMOS observations is that they allow direct 
measurement the polarization fraction as the ratio P/I 
(Figure [2j bottom center). 



The percentage polarization revealed this way varies 
strongly around AB Aur: The disk's southeastern half 
is much less polarized than the opposite side, with 
polarizations of ~ 25 ± 5% and 45 ± 12%, respectively. 
The two bright "shoulders" are seen to be regions 
of maximum polarization (~ 55%), separated by the 
lower polarization (~ 40%) region corresponding to the 
"depleted region" from Opp08. The symmetry axis of 
the overall polarization pattern is 328° ±3°, aligned with 
both the "gap" , and with the inferred rota tion axis of the 
disk HCorder et alj|2005t|p"ietu e t al.ll2005ft . The claimed 
"gap" in polarized intensity is now seen to be a region of 
lower polarization fraction, and not a region of decreased 
total disk-scattered light. There is no significant decrease 
in surface brightness in any of the total intensity images 
at this location (see the white circles in Figure [5]). 

4. A MODEL FOR SCATTERING FROM AB AUR 

The alignment of the polarization pattern with the 
disk's inclined rotational axis provides a convincing 
indication that scattering geometry and not disk sub- 
structure predominantly causes the observed appearance 
in polarized light. As shown in Fig. 1, the polarization 
of dust-scattered light is maximized for scattering angles 
near or just above 90°, and then decreases for scattering 
angles closer to 180°. The higher polarization of the 
northwestern half of the disk indicates that side is tilted 
away from us; the most distant part of the disk, where 
the light is most strongly back-scattered, has lower 
polarization without any need to invoke clearing. 

To demonstrate that this scenario accounts for the 
observed polarization around AB Aur, we produced 
model images u sing the MCFOST Monte Carlo radiative 
transfer code (jPinte et al.l I2006T ) . We concentrated 
on fitting the 2.0 fim polarization by varying dust 
properties, holding most disk parameters (e.g. radii, 
scale height, etc.) fixed bas ed on previous modeling 
(e.g. pTannirku lam et ail 120081 . and references therein) . 
The inclination of AB Aur's disk is not precisely known: 
Using CO kinematics, som e authors have found incli- 
natio ns as high as 33-40° (IPietu et all I2005t ILin et all 
120061 ). but others favor 21° or less (|Corder et all 12005ft . 
Lower inclinations are also supported by the observed 
very low column density of hydroge n and high Ly 
a wind velocity (iRoberee et al.l 12001ft. a nd by near- 
infrared interferometry (lEisner et al J 12004 ). Due to the 
uncertainty, we allowed the inclination to vary as a free 
parameter. 

We evaluated models across a grid in parameter 
space for plausible dust populations (see Table QJ. 
For each model, we computed the SED and 2 fim 
total and polarized images at a range of inclinations. 
Images were convolved with a 0.22" Gaussian to match 
the data's resolution. To find the best fit, we first 
discarded models which did not provide an acceptable 
match to AB Aur's SED, then calculated each model's 
azimuthal polarization profile in an 0.7-1.5" annulus, and 
computed the \ 2 statistic relative to the observed 2.0 /im 
polarization profile. 

The best fit model uses silicate grains with 60% poros- 
ity, maximum grain size 1 fim and a size distribution with 
power law index = —3.5. See Figures [3] and [4] These 
parameters are similar to those inferred for solar system 
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Fig. 2. — Near-IR imaging and polarimetry of AB Aur. Left four panels: our new NICMOS observations. These panels display the 1.1 
and 2.0 fim total intensities (using arcsinh stretches), 2.0 (im polarized intensity (using a log stretch matching that used by Opp08), and 
polarization fraction (shown using a linear stretch indicated by the inset color bar, and with vectors showing the polarization angles). All 
images have been rotated to align the disk's minor axis at PA=328° vertical. Right panels: H band total intensity and polarized intensity 
from previous works, for comparison. 

Our 2 (im polarized intensity observations reproduce the appearance seen by Opp08 very closely, but our P/I image shows that the "depleted 
region" , indicated with a white dashed circle in all panels, is in fact a region of lower polarization fraction, not lower total intensity. None 
of the total intensity images show any indication of an open region at that location in the disk (though we caution that in the H band 
image from Fukagawa et al., the region of interest straddles the radius, shown in black, inside of which they state PSF subtraction artifacts 
rendered their data not photometrically reliable.) 




2 1 -1 -2 2 1 -1 -2 2 1 0-1-2 



Fig. 3. — Our best-fit MCFOST simulated images of AB Aur, for an inclination of 35° and using a population of astrosilicate grains 
with 60% porosity, a max = 1 M m an d dN/da on a -3 / 5 . The three panels show 2 (im total intensity, polarization fraction, and polarized 
intensity, using the same display scales as the corresponding panels in figure [2] See [|4] 
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TABLE 1 

Range of parameters explored in model grid 



90 180 270 

Angle from Disk Minor Axis 



360 



Fig. 4. — Observed and modeled azimuthal profiles for 2.0 /an 
polarization and total intensity, measured in an annulus from 0.7- 
1.5". Azimuthal angles on the horizontal axis are measured relative 
to the minor axis pointing toward us. For both data and models, 
the near side is brightest in total intensity, but the far side is much 
more polarized. The models shown are for the same best fit grain 
population as in Figure 3. For these dust parameters, the double- 
peaked pattern in polarized intensity is best fit at inclinations > 
35°, but total intensity and the overall polarization level is best fit 
at inclinations ~ 22°. We conclude that our model's 60% porous 
silicate spherical grains are only an approximate fit for the dust 
around AB Aur. 



come tary dust grains (jPetrova et all l2000t iShen et all 
l200l . 

We found that inclinations of 22-35° best fit the 
scattered light around AB Aur, consistent with the 
range inferred from CO velocities. With the above 
dust parameters, inclinations of > 35° better reproduce 
the polarization drop from backscattering. Lower 
inclinations < 22° improve the fit to the average 
polarization of the near and far sides, but do not show 
as strong a decrease in polarization on the far side as is 
observed. Lower inclinations also better simultaneously 
fit the total intensity profile as well as polarized intensity 
(See Figured]). 

These slight discrepancies are most likely due to 
oversimplifications in our dust model: MCFOST assumes 
simple Mie scattering from spherical dust grains, an 
imperfect approximation for actual circumstellar dust, 
which is believed to co nsist of irregular fractal aggregates 
(|Dominik et al.l [20071 . and references therein). Realistic 
scattering properties for fractal grains can be calculated 
through more computationally inten sive techniques such 
as the discrete dipole model (e.g. iPetrova et al.l 2000; 



Parameter 



Values 



Inclination 

Dust composition 

Grain porosity 

Grain maximum radius 

Grain size power law 



13—89° , 20 steps evenly spaced in cos i 

astrosilicates a , olivine 1 ", 

ISM silicates+carbon mixture 

0, 0.2, 0.4, 0.6, 0.9 

1, 20, 200 
-2.5, -3.5, -4.5 



Note. - - Fixed model parameters: Stellar T e ff = 9772 K, 
Lt = 471/0, dust mass 10 — 4 Mq, disk Ri nner = 0.2 AU, Router = 
350 AU, scal e height h = 14 AU at R = 100 AU, disk flaring 
h oc i? 1 . Sec Tannirkulam et al. (2008) and references therein. 

a IDraine fcTei lfl98l ) : IWeingartner &: Drainel pOOlT ) 
b IDorschner et al l J1995I) 
c IMathis fc"Wh"i ffcn (1989) 



IDas et alJ200l IShen et al.H2009T) . Ag gregate grains seem 
promising candidates to provide a better fit for AB Aur. 

Specifically, discrete dipole calculations in some cases 
predict polarizations that peak at lower scattering angles 
(90° or below) and decrease more strongly toward 180° 
compared with Mie results (see for instance Figure 
8 of Shen et al. 2009). That shift toward higher 
polarization at lower scattering angles would broaden the 
azimuthal polarization profile and could deepen a dip 
from backscattering, two changes that would improve 
our fit to AB Aur. (In Figure 4, the visible offset 
between the observed profile's polarization peaks, near 
90 and 270°, and the model profiles' peaks, near 135 
and 225°, directly shows the need for dust grains whose 
polarization maximum occurs at lower scattering angle.) 
Improving disk models to use more sophisticated dust 
treatments during radiative transfer is a logical next step. 

5. DISCUSSION: 
5.1. Structure in AB Aurigae's Disk 

The pattern of polarization resolved around AB Aur 
indicates unambiguously that the spatial variation of 
polarized light there is due primarily to the geometry 
of scattering from the inclined disk's surface. These data 
do not support the hypothesis of significant clearing in 
the disk near PA=333°. 

Independent of but simultaneously with this study, 
a numerical investigation of protoplanet shadows in 
disks reached a similar conclusion, finding the observed 
morphology around AB Aur inconsistent with the 
presence of any protoplanet above 0.3 Jupiter masses 
(Jang-Condell and Kuchner, submitted). 

While the previously-claimed gap does not appear to 
be present, there is real structure within AB Aur's disk 
which is visible in scattered light, including polarization. 
Many of the asymmetric spots seen in the 2.0 /jm 
polarization fraction image coincide with the spiral 
arms as seen in the 1.1 /mi image. For instance, 
the northeasternmost region of high polarization, near 
(1.7, 0) in the coordinate system of Figure 2, is precisely 
aligned with the brightest spiral arm. Overplotting or 
blinking these images shows several such alignments. 
Perturbations in the disk's surface, from either localized 
changes in scale height or warps in midplane location, 
could change the scattering geometry at the optical depth 
r = 1 surface to produce the observed minor variations 
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in polarization. For optically thick disks, surface features 
seen in scattered light do n ot necessarily correlate wit h 
conditions at the midplane (|Jang-Condell fc B oss 2007), 
but in the case of AB Aur there is direct evidence 
from CO emission that the b ulk of the gas deviates 
from pure Keplerian rotation (|Lin et alj 120061 ). These 
perturbations and/or warps might contribute to the 
difficulty in firmly establishing the system's inclination. 
But this just raises the inevitable next question: what 
causes those fluctuations? 

Though we find no direct evidence for any disk gap 
due to a protoplanet, there still remains a case for the 
presence of a companion somewhere around AB Aur: 
Both the strong spiral structure seen in dust-scattered 
light and the non-Keplerian dynamics revealed by gas 
emission lines argue for the existence of a planetary-mass 
body perturbing the disk. Other perturbation mecha- 
nisms, such as gravitational instabilities in the disk or a 
stellar-mass companion, seem ruled out observationally 
([Pietu et all [2003: iLin et all 120061 ) . The detection of 
this companion thus awaits future improvements in high- 
contrast imaging to reveal. 

5.2. Interpreting Features in Imaging Polarimetry 

Differential polarimetry, as used by Opp08 and 
others, is a proven technique for observing circumstellar 
dust at high contrast from the ground with adaptive 
optics (AO). Such observations will become increasingly 
com mon with upcoming ex treme AO syst ems, such as 
GPI ([Macintosh et al.l l2006) and SPHERE (|Beuzit et all 
2006). For many disks, these instruments will provide 
high contrast images only in polarized light. 
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We have shown here that care is required when 
interpreting such data. Variation of polarization with 
scattering angle cannot be neglected, even for low 
inclination targets like AB Aur. Yet AO differential 
polarimetry can still yield precise measurements of disk 
structure and dust properties, provided the degeneracy 
between polarization fraction and intensity can be 
broken. This may be possible through multiwavelength 
observatio ns, which are differen tly sensitive to scattering 
geometry ([Watson et al.ll2007t ). or in combination with 
independent constraints on geometry such as from CO 
emission. 

For brighter disks, extreme AO should allow accurate 
measurement of polarization fraction from the ground, 
but for now, HST NICMOS offers the highest polarimet- 
ric precision for measurements of disk-scattered light. 
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